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Muon production at forward rapidity (1.5 < |?7| < 1.8) has been measured by the PHENIX 
experiment over the transverse momentum range 1 < pr < 3 GeV/c in ^/s = 200 GeV p+p coUisions 
at the Relativistic Heavy Ion CoUider. After statistically subtracting contributions from light hadron 
decays an excess remains which is attributed to the semileptonic decays of hadrons carrying heavy 
flavor, i.e. charm quarks or, at high pT, bottom quarks. The resulting muon spectrum from heavy 
flavor decays is compared to PYTHIA and a next-to-leading order perturbative QCD calculation. 
PYTHIA is used to determine the charm quark spectrum that would produce the observed muon 
excess. The corresponding differential cross section for charm quark production at forward rapidity is 
determined to be dcrcc/dj/|y=i.6 = 0.243 ±0.013(stat.) ±0.105(data syst.) (PYTHIA syst.)mb. 

PACS numbers: 13.85.Qk, 13.20.Fc, 13.20.He, 25.75.Dw 
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I. INTRODUCTION 

Measurements of heavy quark production in proton- 
proton {p + p) interactions at collider energies serve as 
important tests for perturbative Quantum ChromoDy- 
namics (pQCD). Bottom production at the Tevatron col- 
lider (V^ = 1.8 and 1.96TeV/c) HQ is reasonably weU 
described by a recent Fixed Order Next-to-Leading Loga- 
rithm (FONLL) calculation 3, 4, 5]. Charm production 
at FNAL, which has only been measured at relatively 
high PT (> 5 GeV/c), is « 50% higher than the FONLL 
prediction Q . However, theoretical and experimental un- 
certainties are large, such that significant disagreement 
between theory and data cannot be claimed. 

Measurements at Brookhaven National Laboratory's 
Relativistic Heavy Ion Collider (RHIC), by both the 
PHENIX and STAR experiments, have provided a wealth 
of information on mid-rapidity open charm produc- 
tion in collisions at ^snn = 130 GcV {p + p) and 
^JJWn = 200 GeV {p + p, d + Au, and Au + Au) 
down to pt ~ 0.5GeV/c. Semileptonic decay of pro- 
duced charm quarks is the primary source of high pT 
leptons after contributions from known (light hadron) 
sources are subtracted. Both PHENIX 0, i, H 113, III 
m m El El 13 and STAR [H El have made sta- 
tistical measurements of charm production via single- 
electron spectra. STAR has also made a direct mea- 
surement of charm production through reconstruction of 
hadronic decay modes of D mesons :17]. In p + p colli- 
sions at y^s^N = 200 GeV PHENIX finds dacc/dy\y=o = 
0.123 ± 0.012(stat.) ± 0.045(syst.) mb U\. STAR finds 
a somewhat higher central value, dacc/dy\y=o — 0.30 ± 
0.04(stat.) ± 0.09(syst.)mb T?!, but the two measure- 
ments are consistent within the stated errors. Both mea- 
surements are noticeably (2-4 x) higher than PYTHIA 
(a leading order pQGD event generator) 0, see ex- 
perimental references for specific parameter sets] and 
FONLL 20]. Again, quantitative disagreement cannot 
be established with current experimental and theoreti- 
cal errors. However, we note that there is some debate 
about whether charm quarks are heavy enough to be re- 
liably treated by pQCD [2ll ]. 

Such measurements also serve as an important base- 
line for charm production in proton-nucleus or deuteron- 
nucleus {p + A or d + A), and nucleus- nucleus (A + B) 

collisions mm mill. In the absence of any nuclear 

effects, charm production (since it is a point-like process) 
is expected to scale with the number of binary nucleon- 
nucleon collisions (Ncoii), which depends on the impact 
parameter of the nuclear collision and can be obtained 
from a Glauber calculation 26]. The degree of scaling 
for any given centrality bin is quantified by the nuclear 
modification factor: 



Rai 



1 ^ dN^^/dy 



Kifi dNPP/dy 



(1) 



Deviations from this scaling {Rab 1) inp + Aoi d + A 
collisions quantify cold nuclear matter eH'ects (such as 
initial state energy loss S IHIHIalEl HI , and shad- 
owing Isl I33. I35I Isl l37j). Any such deviation must 
be understood so that in A + B collisions contributions 
to Rab 7^ 1 from hot nuclear matter effects (such as in- 
medium energy loss 13 §1 and references therein] ) and cold 
nuclear matter effects can be disentangled. In d + Au col- 
lisions both PHENIX and STAR find little or no effect 
of cold nuclear matter on char m p roduction {RdAu ~ 1 
over the measured lepton pT [ifl E3)- This contrasts 
with measurements of open charm in Au -\- Au collisions: 
although the total charm production appears to scale 
with NcoU 0, there is a strong suppression of lepton 
spectra for > 2GeV/c that increases with central- 
ity [ill Iisl ] . Furthermore the elliptical flow of non- 
photonic single electrons, as measured by PHENIX in 
Au + Au collisions [l4L flsL fig] , implies that the charm 
quarks interact strongly with the created medium. 

Finally, since the initial formation of open and closed 
charm are both sensitive to initial gluon densities [3^l40l] . 
open charm production serves as an appropriate nor- 
malization for J /ip production. The production of J /'ip 
mesons is expected to be sensitive to the production of 
a quark gluon plasma (QGP), should one be formed in 
A + B colhsions 0, [4^ El E El . In order to un- 
derstand J /ip production differences in A + B collisions 
compared to p -t- p and p + A collisions it is important 
to take into account any differences in the charm quark 
production in each of the different systems. 

Until now, open charm measurements at RHIC have 
been limited to mid-rapidity. Measurements at forward 
rapidity are interesting for a variety of reasons. First is 
the need to constrain theoretical calculations over a wide 
kinematic range. The importance of this is demonstrated 
by the DO measurement of bottom production at large ra- 
pidity (Vi = 1.8 TeV, PT > 5GeV/c, 2.4 < < 3.2), 
as deduced from the production of high pT muons jj]- 
Significant theoretical improvements resulted from the 
effort to reduce what was, initially, a discrepancy be- 
tween theory and experiment that increased with increas- 
ing rapidity @. Second, significant cold nuclear effects 
have been seen in RHIC collisions at forward rapidity. 
PHENIX 113, BRAHMS ^ ^], and STAR ^ have 
all measured light hadron production in d + Au collisions 
at forward rapidity and have found significant deviations 
from RdAu = 1- It will be interesting to see whether 
charm production follows a similar pattern. Finally, open 
charm production at forward rapidity needs to be under- 
stood to fully interpret PHENIX J/ip measurements at 
forward rapidity 0, (H El El ■ 

In this paper we report on the measurement of muon 
production at forward rapidity (1.5 < \ri\ < 1.8), in the 
range 1 < pr < 3 GcV/c, in ^ = 200 GcV p + p col- 
lisions by the PHENIX experiment. The upper limit 
of the Pt range is determined by available statistics. 
The vertex-independent muon yield is statistically ex- 
tracted by calculating and subtracting contributions from 
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light mesons (tt's and ivT's) which decay into a muon, 
and hadrons which penetrate through the muon arm ab- 
sorber material. In the absence of new physics, and 
in the px range measured in this analysis, such muons 
come dominantly from the decay of hadrons containing 
a charm quark (with small contributions from decays of 
hadrons containing a bottom quark and decays of light- 
vector mesons) . PYTHIA is used to determine the charm 
quark spectrum that would produce the observed vertex- 
independent muon spectrum, and from this we obtain 
the differential cross section of charm quark production 
at forward rapidity. 

The remainder of this paper is organized as follows: In 
Section we describe the PHENIX experimental appa- 
ratus, with special emphasis on the muon arm detectors. 
In Section ITTll we describe the methodology used to ex- 
tract the vertex-independent muon signal. This section 
includes details on the run, event and track selection cri- 
teria; values obtained for contributions to the muon yield 
from abundant light hadrons, which are subtracted to 
obtain the vertex-independent muon yield; and details 
on the systematic error analysis. In Section IIVI we ex- 
tract the differential cross section for charm production 
at y — 1.6, integrated over pr- Finally, in Section Ivl 
we compare to other measurements, draw conclusions, 
and discuss the prospects for such measurements with 
improved data sets currently under analysis. 



II. THE PHENIX EXPERIMENT 

The PHENIX experiment , shown in Figure ^ is 
a large multipurpose set of detectors optimized for mea- 
suring relatively rare electromagnetic probes (photons, 
muons, and electrons) of the spin structure of the proton 
and of the hot dense matter created in ultrarelativistic 
heavy ion collisions. The data acquisition system and 
multilevel triggers are designed to handle the very dif- 
ferent challenges presented hy p + p collisions (relatively 
small events at very high rates) and Au + Au collisions 
(very large events at relatively low rates) with little or no 
deadtime ,5,5, 56|. Event characterization devices, such 
as the Beam-Beam Counters used in this analysis, 
provide information on the vertex position, start time, 
and centrality of the collision. The two muon arms cover 
1.2 < I77I < 2.4 in pseudorapidity and 54> = 27r in az- 
imuth. The two central arms, which each cover |?7| < 0.35 
and (50 = 7r/2, are not used in this analysis. 

The Beam-Beam Counters (BBCs) 57] each consist 
of 64 quartz radiators instrumented with mesh dynode 
PMTs and arranged in a cylinder coaxial with the beam. 
The BBCs are placed on either side of the collision ver- 
tex and cover 3.0 < \r]\ < 3.9. Each channel has a dy- 
namic range extending to 30MIPs. The BBCs measure 
the arrival times of particles on both sides of the colli- 
sion vertex, ts and t^. From the average of these times 
we determine the event start time. From their difference 
we obtain the position of the vertex along the beam di- 



rection, Zytx- The BBCs also provide the minimum bias 
interaction trigger, which requires that there be at least 
one hit in each BBC and that \zvtx\ < 38cm. 

The muon arms |58j | are coaxial with the beam on op- 
posite sides of the collision vertex. By convention the 
arm on the South (North) end of the interaction region 
is assigned negative (positive) z coordinates and rapid- 
ity. For the 2001/2 run period, in which the data for 
this paper were collected, only the South muon arm was 
operational. Each muon arm is comprised of a Muon 
Tracker (MuTR) and a Muon Identifier (MuID). The 
MuTR makes an accurate measurement of particle mo- 
menta. The MuID allows coarse resolution track recon- 
struction through a significant amount of steel absorber. 
Together the muon arm detectors provide significant pion 
rejection (> 250 : 1, increasing with decreasing momen- 
tum) through a momentum/penetration-depth match. 

Before ever reaching the MuTR detectors a particle 
must pass through the pre-MuTR absorber: 20 cm of 
copper (the nosecone) plus 60 cm of iron (part of the 
MuTR magnet). The nominal nuclear interaction lengths 
of iron and copper are Xf'^ — 16.7cm and Xf'^ — 15.3cm 
(although this varies with particle species and energy, 
see Section 1111 !*'() ■ Therefore the pre- MuTR absorber 
presents a total thickness of 4.9A// cos 9, where 6 is the 
polar angle of a particle's trajectory. This absorber 
greatly reduces the MuTR occupancy and provides the 
first level of pion rejection. 

Each MuTR arm consists of three stations of cathode 
strip chambers installed in an eight-sided conical mag- 
net m. The radial magnetic field ( / B • dl = 0.72 T-m 
at 15 degrees, B{9) « 5(15°) tan(0)/ tan(15°)) bends 
particles in the azimuthal direction. Each station occu- 
pies a plane perpendicular to the beam axis and consists 
of multiple ionization gaps (3 gaps for the two stations 
closest to the collision vertex, 2 gaps for the last station) 
which have their charge imaged on two cathode strip 
planes oriented with a small stereo angle to provide two- 
dimensional information. An ionizing particle typically 
fires three adjacent strips in each orientation. A fit to the 
charge distribution on adjacent strips provides a position 
measurement with a resolution of cr ~ 100 /im in the bend 
direction. The MuTR achieves a momentum resolution 
of Cp/p « 5% over the analyzed kinematic range. The 
resolution is approximately independent of momentum 
due to the significant contribution from energy loss fluc- 
tuations in the pre-MuTR absorber, which falls as 1/p, 
and which counters the more familiar linear momentum 
dependence seen for particles tracked through a "thin" 
detector. 

Each MuID arm consists of five steel absorber plates 
interleaved with larocci tubes (operated in proportional 
mode) and specialized shielding to reduce backgrounds 
not originating from the collision vertex. Gaps are 
labeled 0-4 proceeding downstream from the collision 
point. 

The larocci tubes, which are between 2.5 and 5 m in 
length, have eight 1 cm^ square cells, each consisting of a 
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FIG. 1: (Color online) PHENIX experimental layout during the 2001/2 run period. 



three-sided ground electrode and an anode wire, mounted 
inside a PVC gas enclosure. A readout channel ("two- 
pack") is formed by wire-ORing the 16 anode wires of 
two tubes which are mounted in planes perpendicular to 
the beam axis and staggered by half of a cell width (0.5 
cm). This provides redundancy, eliminates geometric in- 
efficiency due to the cell walls, and reduces the maxi- 
mum drift time for charge collection. Digital readout of 
the two-pack signals provides a coarse one-dimensional 
hit position {a = 9cm/\/T2 = 2.6 cm). The tubes in 
each gap are mounted in six individual panels, each of 
which contains two layers of two-packs (horizontally and 
vertically oriented), thus providing two-dimensional in- 
formation. 

The first MuID absorber plate (thickness = 20 cm - 
South; 30 cm - North) also serves as the return yoke 
of the MuTR magnet. Successive plates (identical for 
the two arms) are 10, 10, 20 and 20 cm thick, thus to- 
taling 4.8A//cos6' (5.4A//cos6') for the South (North) 
arm. Due to ionization energy loss a particle must have a 
momentum at the vertex which exceeds 2.31 cos0 GeV/c 
(2.45 cos6' GeV/c) to penetrate to the most downstream 
MuID gap of the South (North) arm. 

Steel plates surrounding the beam pipe guard against 
backgrounds caused by low-angle beam-beam collision 
products which scrape the beam pipe near the MuID 
2;-location (7-9 m) or shine off the RHIC DX magnets 
immediately downstream of each MuID arm. Steel 
blocks in the RHIC tunnels guard against penetrating 
radiation generated by the incoming beams scraping 
against beamline components (primarily the final focus- 



ing quadrupoles). 

The MuID also contributes information to the first- 
level trigger decision. For the 2001/2 run, during which 
the data for this analysis were collected, a relatively 
coarse trigger was implemented using LeCroy 2372 Mem- 
ory Lookup Units (MLUs). Each gap was divided into 
quadrants with a horizontal and vertical split going 
through the beam axis. Signals from tubes in an individ- 
ual gap/orientation (layer) and quadrant were logically 
combined. Only gaps 0,2,3 and 4 were used in the trig- 
ger due to the 16-bit input limitation of the MLUs. The 
penetration depth required for the trigger to fire was pro- 
grammable. The MuID-lDeep trigger fired if more than 
6 out of 8 layers in a particular quadrant were hit (indi- 
cating the possibility that the event contained a particle 
penetrating to MuID gap 4). The MuID- 1 Shallow trig- 
ger fired if 3 of the 4 most shallow layers (horizontal and 
vertical layers in gaps and 2) were hit for a particular 
quadrant. 



III. METHOD FOR EXTRACTION OF MUONS 
FROM CHARM DECAY 

Inclusive muon candidates, Nj, are those particles 
which are successfully reconstructed to the last MuID gap 
(gap 4). These consist of four components: 1) "free-decay 
muons" , Nu , which result from the decay of light hadrons 
(tt and K mesons) before reaching the pre-MuTR ab- 
sorber, 2) "punchthrough hadrons" , Np, which penetrate 
the entire detector and are thus misidentified as muons 
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FIG. 2: (Color online) Schematic depiction of the relative flux of different components of the inclusive muon candidate yield 
as a function of flightpath into the muon arm absorber (the event vertex is at z-utx = 0). See text for details. 



3) "background tracks" , Nb , which in p + p collisions 
are dominated by hadrons which decay into a muon after 
reaching the MuTR, and 4) "vertex-independent muons" , 
N^, which are primarily due to the decay of heavy flavor 
mesons. 

Figure El shows a schematic depiction of the relative 
yield per event of these different contributions as a func- 
tion of flightpath into the muon arms, as described below. 

The number of hadrons is large and essentially in- 
dependent of flightpath until the first absorber layer is 
reached. In each absorber layer these hadrons undergo 
strong interactions with a probability 1 — exp(— L/A), 
where L is the length of absorber material traversed, 
and A is the species and pT-dependent nuclear interac- 
tion length determined in Section IIII Fl Most of these 
interacting hadrons are effectively eliminated as pos- 
sible muon candidates. However, a small fraction of 
hadrons penetrate the entire absorber thickness. These 
punchthrough hadrons are indistinguishable from muons. 

The decay lengths for tt's (cr = 780 cm) and if's 
(cr — 371 cm) are long compared to the flightpath from 
the vertex to the absorber. Therefore, the fraction of de- 
cay muons from these sources is relatively small, but in- 
creases linearly with the flightpath until the first absorber 
layer is reached. A hadron which decays prior to the pre- 
MuTR absorber into a muon that is inside the detector 
acceptance is indistinguishable from a muon originating 
at the vertex. After the first absorber layer the number 
of free-decay muons remains constant by definition. 



Hadrons which decay in the MuTR volume are a rel- 
atively small contribution since most are absorbed prior 
to reaching the MuTR, the Lorentz-dilated decay lengths 
are long compared to the length of the MuTR volume 
(South « 280 cm, North 420 cm), and a particle which 
decays in the MuTR is less likely to be reconstructed. 
Such tracks are partially accounted for in the calculation 
of punchthrough hadrons (see Section Fill F|l and the re- 
maining fraction falls under the category of background 
tracks (Section IIII G|l . This small contribution is not 
shown. 

Without a high-resolution vertex detector muons from 
various sources (the decay of open heavy flavor hadrons, 
the decay of quarkonia, the decay of light vector mesons, 
and Drell-Yan production) originate indistinguishably 
close to the collision vertex. Thus their yield is inde- 
pendent of the flightpath and independent of the vertex 
position. Since inclusive muon candidates, by deflnition, 
penetrate to MuID gap 4, we measure the combined yield 
at 2 w 870 cm. 

Figure 13 shows a sample distribution of the inclusive 
muon candidate yield as a function of collision vertex 
{zytx), and its decomposition into the four different con- 
tributions. The yield of free-decay muons is seen to have 
a linear dependence that is set to at z^tx — Zabs — ^d- 
Here Zabs = —40 cm is the upstream face of the pre- 
MuTR absorber (indicated by the thick solid line), and 
A^) is the effective absorption length, beyond which there 
are no free-decay muons. A^i was found to be nearly iden- 
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FIG. 3: Sample z^tx distribution of different components 
of the inclusive muon candidate yield (measured data for 
both charge signs over the range 1.0 < pr < 1.2GeV/c). 
Crosses show inclusive muon candidates, filled circles show 
free-decay muons, open circles show punchthrough hadrons, 
open squares show background tracks, and open diamonds 
show the sum of these three hadronic sources. The vertex- 
independent muon yield is obtained from the difference be- 
tween the yield of inclusive muon candidates and the yield of 
light hadronic sources. 



tical to the species and p-r-dependent nuclear interaction 
lengths determined in Section IIII Fl Muons originating 
from meson decays downstream of this location have no 
Zytx dependence. The fraction not accounted for in the 
calculation of the punchthrough hadron yield forms the 
small contribution from background tracks. The yield 



of punchthrough hadrons and vertex-independent muons 
also have no Zytx dependence. Note that the ratio of 
different contributions to the inclusive muon candidate 
spectrum is pt dependent. 

In order to extract the cross section for charm pro- 
duction we first need to determine the yield of vertex- 
independent muons, N^{pt), the amount beyond that 
due to light hadrons and fake backgrounds. As described 
in Section IlII Al we select good runs, events, and tracks, 
and restrict our acceptance to regions where the detec- 
tor performance was optimal, and the acceptance vs. z^tx 
was nearly constant. Next, as described in Sections IIII Bl 
and IIII CI we obtain the yield of inclusive muon candi- 
dates vs. Pt and z^tx, corrected for acceptance and ef- 
ficiency: Ni{pT,Zytx)- In Section [ill Dl we describe a 
data-driven hadron generator. This generator is used in 
Section flll El in which we describe how the vertex depen- 
dence of the inclusive muon candidate yield allows us to 
obtain the yield of muons from light-meson decay before 
the pre-MuTR absorber, similarly corrected and binned: 
NoiPT, 



• Zytx)- This generator is also used in Section llll Fl 
in which we describe how we use hadrons which stop in 
MuID gap 3 (the penultimate gap), together with simu- 
lations of hadron penetration in the MuID absorber, to 
obtain the yield of punchthrough hadrons in MuID gap 
4: Np{pT, Zytx)- The yield of fake tracks, Nb{pt, Zvtx), 
determined from simulations described in Section IIII Gl 
is found to be small. 

The yield of vertex-independent muons is determined 
by subtracting the contributions from light hadrons and 
fake backgrounds and averaging over Zytx bins: 



Ni{pT, zitx) - Nd{pt, zlt^) - Np{pT, - NsipT, 



(2) 



where (P /2'KpTdpTdri is implicit in all terms of the equa- 
tion. 

We convert this into a cross section via 

2TTpTdpTdy esse' '^''^PTdprdr] 

Here (J^g^Q is the cross section of the BBC trigger for 
p + p interactions and e^'g^'' is the efficiency of the BBC 
trigger for events in which a charm quark is created and 
decays into a muon. Substituting ?7 — > y introduces neg- 
ligible error due to the small mass of the muon, the only 
component left after the subtraction. As described in 
Section IIII II systematic errors are determined for each 
component and combined into a term that applies to the 
overall normalization and a term that applies to the pT 
dependence of the spectrum. 



We use PYTHIA to derive the pT-dependent differ- 
ential cross section for the production of charm quarks 
responsible for the vertex-independent muon yield. This 
procedure is very similar to that in references 0, IE IS 
llOl IllL [12. , and is described in detail, along with the 
associated systematic error analysis, in Section HVI 



A. Data Reduction 

1. Data Sets and Triggering 

Runs were selected for this analysis based on stable 
detector operation using the same criteria as an earlier 
analysis of J/ip production |5ll |. Further runs were elim- 
inated due to the presence of large beam-related back- 
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grounds entering the back of the detector. 

We select only those events in the vertex range — 20 < 
Zvtx < 30 cm. This minimizes the Zytx dependence of the 
detector acceptance and allows us to treat the amount 
of absorber material as a simple function of polar angle, 
ignoring complications in the pre-MuTR absorber near 
the beampipe. 

The decision to collect an event was made by the Lo- 
cal Level-1 Trigger [s^ within 4 /is of the collision. Input 
to the trigger decision was given by the BBC (collision 
with a valid event vertex) and the MuID (reconstructed 
penetrating track). Each trigger could be independently 
scaled down, so that it occupied a predetermined frac- 
tion of the available bandwidth, by selecting every N^^ 
instance, where Ni is the scaledown factor for the i*'' 
trigger. Three different data sets were selected from the 
good runs for different aspects of the data analysis: 

• BBC: To extract Njj we need to measure the z^^j. 
dependence of Nj. For this we need the unbi- 
ased collision vertex distribution, which we obtain 
from a set of events collected with the BBC trigger: 
Nbbcn > 1 && Nbbcs > 1 && \zvtx\ < 38cm, 
where Nbbcn and Nbbcs are the number of hits 
in the North and South BBC respectively, ct^bbc 
was found to be 21.8 ± 2.1 mb using a van der Meer 
scan |60| . There were 1.72 X 10^ BBC triggered 
events passing our vertex selection criteria in this 
data set, corresponding to a sampled luminosity of 
J Ldt = 0.79 nb-i. 

• MuID- Weep && BBC {MID): In order to extract 
Ni, Nd and Nb we used events selected with the 
MID trigger, which enriched the sample of events 
with tracks penetrating to MuID gap 4. For the 
MID and BBC data sets we used identical run se- 
lection criteria. The total number of sampled BBC 
triggers for this data set was 5.77 x 10^, correspond- 
ing to a sampled luminosity of J Ldt = 26.5 nb~^. 

• MuID-lShallow kk BBC (MIS): In order to ex- 
tract Np we need a data set which provides an unbi- 
ased measurement of the number of particles which 
penetrate to MuID gap 3. Since the MID trig- 
ger required tracks to penetrate to MuID gap 4 it 
could not be used. Instead we used the Mi 5* trigger, 
which only used information from MuID gaps 0-2. 
We used a subset of runs for which the scaledown 
factor for this trigger was only 10, corresponding 
to a sampled luminosity of J Ldt — 1.72 nb^^. 



2. Track Selection 

The Muon arm reconstruction algorithm starts by 
finding "roads" (collections of hits in the MuID which 
form straight, two-dimensional lines) and then combin- 
ing them with hits in the MuTR to form "tracks". We 
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FIG. 4: The scaled angular deflection is the difi'erence in a par- 
ticle's polar angle caused by passage through the pre-MuTR 
absorber scaled by the particle's momentum, p69. For muons 
(and hadrons not undergoing a strong interaction in the pre- 
MuTR absorber) one expects the distribution of this quantity 
to be well described by the standard multiple scattering for- 
mula. The different panels show p59 for different pt bins with 
fits (normalization only) to the expected distribution. 



apply strict cuts on both road and track parameters in 
order to reduce backgrounds, see Tabled 

The resulting purity of the selected tracks is demon- 
strated in Figure ^ This figure shows pS9, the angu- 
lar deviation through the pre-MuTR absorber, scaled by 
the particle momentum to give a quantity which should 
be momentum independent, for different pt bins. As 
shown in Figure El S9 is the angular difference between 
the reconstructed particle trajectory at the collision ver- 
tex {x = 0, y = 0,z = Zytx) and at MuTR station 1. A 
GEANT Hi] simulation of the PHENIX detector showed 
that tracks which do not suffer a strong interaction 
in the pre-MuTR absorber undergo angular deviations 
consistent with expectations based on standard multi- 
ple scattering: aso oc ^Jx/Xq/p. The curves in each 
panel are fits to Cp59 eyi\){—{p59)'^ /2{p<Jse)^), in which 
the normalization constant (C) is allowed to float, and 
pase = 130rad-MeV/c is given by GEANT and is consis- 
tent with a simple estimate based on the radiation length 
of the pre-MuTR absorber and the standard multiple 
scattering formula [Hilil (x/Xq « 48 pOlZle ~ 



9 



TABLE I: Road and track cuts. Here Dp is the penetration 
depth, defined to be the most downstream MuID gap with at 
least one hit (from the horizontal or vertical layer) associated 
with the track; z is the coordinate along the beam; x and y are 
transverse to each other and to the beam axis; the vertex cut 
refers to the transverse position of the MuID road projected 
to the xy plane at 2 = 0; and the slope cut refers to the 
direction cosine of the road in each transverse direction. 

# Associated MuID hits, NmuID > 6 

(out of a possible 2Dp) 
Road cuts Vertex cut, y^x^ + y^ < 100 cm @ z = 

Slope cut, 0t)2 + (f )2 > 0.25 
> 1 associated hit in MuID gap 4 
Track cuts Track fit quality, x'^/dof < 10 

# Associated MuTR hits, Nmutb. > 12 

(out of a possible 16) 
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FIG. 5: The angular deflection, 59 is the angular difference 
between the reconstructed particle trajectory at the collision 
vertex and at MuTR station 1. The momentum used to scale 
56 is the average of the momentum reconstructed inside the 
MuTR magnet (pstai) and the momentum extrapolated to 
the vertex (pvtx)- 

(V2)(13.6MeV/c)(^/48)rad 133rad-MeV/c). The in- 
tegral beyond Spase is ~ 5% and is largely due to hadrons 
which have a strong interaction in the pre-MuTR ab- 
sorber and are still reconstructed as a muon candidate. 
Such tracks are accounted for in the calculation of the 
punchthrough hadron yield, as described below. 



3. Acceptance Restriction 

We further restricted the acceptance of muon candi- 
dates for this analysis in two ways: 

1. We required tracks to pass through 6/(p regions in 
which the Monte Carlo detector response strictly 
agreed with the measured response. This was es- 
tablished by agreement between the number of data 



hits and Monte Carlo hits assigned to tracks in each 
9/(1) region of the detector. 

2. We required tracks to lie within a pseudorapidity 
range, 1.5 < |ry| < 1.8, a region over which the ac- 
ceptance depends only weakly on the collision z^tx 
location. 



B. Acceptance and Efficiency 

We factorized the acceptance and efficiency for tracks 
penetrating to a particular MuID gap, i, into four com- 
ponents: 

1. e!^j.j.: the acceptance of a perfectly working de- 
tector with the acceptance restrictions described 
above. This quantity (w 50%) is normalized to 
2^61^ [Srj — 0.3) and accounts for non-sensitive 
structural members in between the cathode strip 
chambers and chamber regions removed from con- 
sideration for the entirety of this analysis. 

2. e*g(.: the efficiency of reconstructing a track that 
fell within the accepted region. This quantity is 
somewhat low (64%) due to detector problems in 
this first physics run that have been subsequently 
resolved. 

3. eJisor: the efficiency for reconstructed tracks to pass 
the cuts listed in Tabled 

4. eJiig: the efficiency of the MuID trigger to fire in 
events with selected tracks. 

e^cc: Srec: ^'^'^ ^user wcrc evaluated with a GEANT 
simulation using single muons thrown with a realistic pt 
spectrum into the muon arms. The applied detector re- 
sponse incorporated measured detector performance. Re- 
ductions in efficiency due to occupancy are negligible in 
p-\-p collisions. Run-to-run variations were ignored since 
we selected runs in which the detector performance was 
similar and stable. Efficiency values for tracks penetrat- 
ing to MuID gap 4 were parameterized in terms of z^tx 
and Pt and are listed in Table ^ There are minor dif- 
ferences in these parameterizatioirs for particles with dif- 
ferent charge sign. 

We also determined the efficiencies for tracks which 
only penetrate to MuID gap 3, since these are needed to 
obtain the yield of punchthrough hadrons. These were 
found to scale from the efficiencies for tracks penetrat- 
ing to MuID gap 4: el^^e^^^^el^^, = el^^y^ x el^^e^^ei^^,, 
where £3^.^10 = 0.66. Egcaie ^^^^ than one because the 
MuID and the road reconstruction algorithm are opti- 
mized for deeply penetrating particles. Particles which 
do not penetrate to the last gap have poorer resolution 
matching to MuTR tracks (due to reduced lever arm and 
smaller number of hits) and are also more susceptible to 
MuID inefficiencies. 
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TABLE II: Trigger, acceptance, track reconstruction and 
track selection efficiencies. Systematic errors for these quan- 
tities are given in Tables IVI and rVITI 



Quantity 



Value 



0.51 X (1 - 114exp(-5.9pT)) x (1 + 0.0015z^t^ 
0.50 X (1 - 531 cxp(-7.5pt)) x (1 + 0.0013z„ti 
0.64 

0.74 X (l-0.00192„t:,) 
0.74 X (1 - 0.0009z„te) 

0.66 

0.86 

0.97 

0.75 



^acc 
^rcc 
user 



p3 

^scalc 



'^BBC 



Trigger efSciencies, e\^-^^ and Eji-ig' ^''^ ^^^o listed in 
Table Ull These were evaluated using the BBC data set, 
which did not require the MuID trigger to fire. 



-trig 



{N^\M1D) X Smid 
{Ni\BBC) X Sbbc' 



(4) 



where Ni\MlD is the number of selected tracks penetrat- 
ing to MuID gap 4 for events in which the MID trigger 
fired, Smid is the scaledown factor applied to the MID 
trigger, and similarly for MID BBC . e^rig ^^^'^ 
evaluated according to Equation 0] but with A^4 N3, 
and MID MIS. 

Since both the MID and MIS triggers required the 
BBC trigger in coincidence with the MuID trigger we 
must also account for the BBC trigger efficiency for 
events in which a reconstructed muon is created via 



charm quark decay: s'bbc' ■ The BBC efficiency was 
evaluated for events in which a J/ijj was produced in 
the muon arm acceptance using PYTHIA-f GEANT sim- 
ulations 

The BBC efficiency was also evaluated for 
events in which tt^'s were produced in the central arm 
acceptance |63| using data triggered without a BBC re- 
quirement. The BBC efficiency under both conditions 
was found to have a similar value that we therefore adopt: 

^BBC — 

Systematic errors for all acceptance and efficiency cor- 
rections are discussed in Section IIII II and listed in Ta- 
bles |V| and EIII 



C. Inclusive Muon Candidates 

We first form two sets of collision vertex {z^tx) his- 
tograms with 10 cm bins: one histogram for all inter- 
actions selected with the BBC trigger, and a series of 
histograms for interactions selected with the MID trig- 
ger and having a good muon candidate within a pt bin 
(1 < PT < 3GeV/c, 200MeV/c bins). The muon can- 
didate histograms are formed separately for each charge 
sign. Entries into each histogram are scaled by the ap- 
propriate trigger scaledown factor. The muon candidate 



histograms are divided by the minimum bias histogram 
to give Ni{pT, Zytx), as shown in Figure El Systematic 
errors shown in this figure are discussed in Section Fill II 
and hsted in Table [V] 



D. Hadron Generator 

In order to determine the contributions to the inclu- 
sive muon yield from free-decay muons f Section IIII E|l . 
and punchthrough hadrons (Section IIIIF|) we make use 
of a data-driven hadron generator. The input for this 
generator is obtained from PHENIX measurements in 
y/s — 200 GeV p + p collisions at y = ,6^ using 
the following procedure: 

1. 7r+ and tt^ spectra at y = (0 < pt < 3.5GeV/c) 
are fit to a power law. We assume factorization in y 
and Pt and scale the spectra fit at y = according 
to: 



^;=obT)exp(-(1.65V2a2)), 

with ffy = 2.5. This factorization is observed both 
in PYTHIA and in BRAHMS ^ data measured 
at y = 1 and y = 2.2. 

2. We use a similar procedure to obtain the charged 
kaon yield at y = 1.65, but we need to first ex- 
trapolate the yield at y = beyond the current 
measurement limit {pr < 2GeV/c). 

We start by forming the isospin averaged K/t: ratio 
vs. Pt &t y = 0. For pT <2 GeV/c we use charged 
particles, [K^ +K^)/{ti^ +it~) [i^- We use neu- 
tral particles, i^°/7r°, for 2 < pt < 6.5GeV/c [H. 
We then fit this combined ratio to the form /{pt) = 
A{1 — B exp{—CpT))- Next, we normalize this 
function separately to the K~^/tt~^ and K~ /tt~ ra- 
tios for PT < 2GeV/c. Finally, we multiply by 
the corresponding charged pion spectrum to ob- 
tain N^^q{pt), our parameterization of the mid- 
rapidity charged kaon pt spectra extending out to 
3.5GeV/c. 

As with pions, we need to extrapolate this param- 
eterization of the yield at y = to obtain the yield 
at y = 1.65. One possibility is to assume boost in- 
variance of the K/n ratio. However, PYTHIA gives 
a slightly narrower rapidity distribution for kaons 
than for pions, resulting in a kaon yield at y = 1.65 
that is only 85% of that predicted with the boost 
invariance assumption. We split the difference be- 
tween these two assumptions: 
^± 

''y=1.65y 

where, again, ay = 2.5. 



iV5, .,{pt) - 92.5%A^^C'o(PT)exp(-(1.65V2a2)), 



3. The p and p spectra are assumed to have the 
same shape as the pion spectra with normaliza- 
tion factors set to the measured values at y = 0, 
prp = 3GeV/c (0.4 for p/n+, 0.24 for p/tt") [6|. 
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FIG. 6: Yield of (left) positively and (right) negatively charged inclusive muon candidates vs. z^tx for different pr bins. Fits 
shown use the functional form a + bzvtx to extract the contribution from hadron decay, as discussed in Section IIII El Error 
bars show statistical errors. Shaded bands show systematic errors, as discussed in Section Hini and listed in Table IVI 



The exact form used for the p, p spectra is unimpor- 
tant. They obviously do not contribute to the yield 
of decay muons and their contribution to the yield 
of punchthrough hadrons is greatly suppressed due 
to their relatively short nuclear interaction length. 

Systematic errors associated with this hadron genera- 
tor are discussed in Section ITlI II and listed in Table EH 



E. Free-Decay Muons 

In Figure El one can clearly see the linear dependence 
in the yield of inclusive muon candidates vs. Zytx at low 
transverse momentum (pT < 2GeV/c). This depen- 
dence is due to muons from the decay of abundant light 
hadrons (tt's and K^s) prior to the first absorber material 
at Zabs — —40 cm. We fit these histograms with the func- 
tion a-\-bzytx- After multiplying by dz/dlfp = cos(< 9 > 
) — 0.947 the slope, b, and its fit error give, respectively, 
the yield per unit length of decay flightpath of muons from 
hadron decay, dN£){pT) / dlfp, and the statistical error on 
this quantity. Results are shown in Figure [7| Systematic 
errors shown in this figure are discussed in Section IIII II 
and hsted in Table |Vl 

This procedure does not provide a quantitative mea- 
sure of the decay muon spectrum above pt 2GeV/c, 



even though a substantial fraction of the inclusive muons 
are decay muons up to significantly higher px- This is 
due to the fact that at high px the decay slopes decrease 
(Lorentz time dilation) as do the statistics, both of which 
make it more difficult to quantify the decay component 
directly. In order to extend our estimate of decay muons 
to higher px we use our hadron generator, described in 
Section lillDI We simulate the decay of hadrons into the 
muon arms and obtain predicted px spectra (per unit 
length) of muons from hadron decay separately for each 
charge sign. We then normalize these predicted spec- 
tra to the measured spectra. The normalized predicted 
spectra are shown as the dashed lines in Figure \7\ The 
predicted spectral shape agrees with the data where we 
have a statistically significant measurement. The abso- 
lute normalization of the prediction is within 7% of the 
measured value, easily consistent within errors. 

We obtain N]j{px, z.^tx) from the product of 
dNx)ipx)/dlfp and the average value of the decay flight- 
path, Ifp = Xd + \zvtx - Zabs\/ cos{6), for each Zytx bin. 



F. Punchthrough Hadrons 

A hadron penetrating to MuID gap 4 is impossible 
to distinguish from a muon. However, we can cleanly 
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identify hadrons in shallower gaps and then extrapolate 
their yield to obtain the yield of punchthrough hadrons 
in MuID gap 4. 

Figure |S1 shows the longitudinal momentum (p^ , the 
momentum projected onto the beam axis) distribution 
of particles that stop in MuID gap 3. The sharp peak at 
Pz ~ 2.2GeV/c corresponds to charged particles which 
stopped because they ranged out in the absorber plate 
between gaps 3 and 4 (this includes both muons and also 
hadrons which only suffered ionization energy loss.) The 
width of this peak is due to the 20 cm (11.4 Xp) absorber 
thickness between MuID gaps 3 and 4, and energy-loss 
fluctuations in all the preceding absorber layers. Parti- 
cles at momenta beyond the peak {p^ > 3GeV/c) form 
a relatively pure sample of hadrons, with only a small 
contamination due to inefficiencies in MuID gap 4 and 
particles with mis-reconstructed momentum values. Af- 
ter correcting for acceptance and efficiency we use these 
particles to obtain the pt spectrum for the "gap 3 exclu- 
sive yield," as shown in Figure^ Note: we use data from 
the MIS trigger sample since the MID sample required a 
hit in MuID gap 4, which would bias this measurement. 

In order to extrapolate this measured spectrum for 



hadrons stopping in MuID gap 3 to the spectrum of 
punchthrough hadrons which penetrate to MuID gap 4 
we start by assuming exponential absorption of hadrons 
entering the muon arm absorber material. With this as- 
sumption we obtain an expression for the "gap 3 inclusive 
yield", those hadrons that reach at least MuID gap 3: 

N^.iPT, e) = Nl,,{pT. 9) expi^L,ie)/yipT)), (5) 

where i indicates the contributing hadron species 
{n^ , ,p,p), Nl^^{pT,0) is the yield at the vertex of 
the i*'' species, ^3(6') is the amount of absorber material 
traversed to reach MuID gap 3, and A'(pt, S) is the pT- 
dependent nuclear interaction length of the i*'' species. 
We can write a similar expression for the punchthrough 
hadron yield: 

N'pipT, 9) = K^^ipT, 9) eM-Lm/\\pT)). (6) 

where L^^ifi) is the amount of absorber material traversed 
to reach MuID gap 4. 

By taking the difference between these two equations 
we obtain an expression for the gap 3 exclusive yield: 



iVr*°^(PT, 9) = Nf^ipT, 9) - iV>(pT, 9) 

= Nl,,{pT.9)eM~Lmi><\PT)) 
x(1 - exp((L3(^?) - Um/>^\VT))) (7) 

In our measurement we cannot identify the species comprising the gap 3 exclusive yield, but we do know their 
charge sign. As a result, Equation [3 can be rewritten as two equations with six unknowns for each pT bin: 

Np^'^^iPT, 9) = N+{pT, 9) - N+{pT, 9) 

- ^ Nl,,ApT.9)eM-Lm/>^'iPT)) 

X (1 - exp((L3(0) - U{9))I\\pt))), (8) 



N-'''°^{pT.i 



x{l-eM{Lm-Lm)/\\pT))). 



N^{pT,9)-Np{pT,9) 

K,^{pT,9)eM~Lsi9)/X\pT)) 



i — TT .K .p 



(9) 



Based on measured cross sections for various 
species [6^ , we chose to reduce the number of unknowns 
with the following assumption: 

Xk+ — Xlo7ig, 
= Ajr+ = A^r- = A/^- = Xshorti 

Xp = 0. 

We further assume that Xshort and Xiong have the form 
a + h{pT[GeY/c] - 1). 



We effectively smoothed the gap 3 exclusive yield for 
each sign by fitting the measured values to a power law. 
Using N^f^{pT,9) from the hadron generator (normal- 
ized to the free-decay muon spectrum, as described in 
Section lill E|l and known values for ^3^4(6*), we fit Equa- 
tions |H1 and |51 to the smoothed gap 3 exclusive yield for 
each sign to obtain: 

Xshort = 19.0 + 2.2(pT[GeV/c] - l)cm, and 
Xiong = 25.9 + 4.4(pT[GeV/c] - l)cm. 
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FIG. 7: Yield per unit length of (left) positively and (right) negatively charged free-decay muons. Points are the measured 
values determined by linear fits to the inclusive muon candidate yield (Figure |SJ. Error bars indicate statistical errors for those 
fits, pt bins with a negative (non-physical) slope in those fits are shown with a line at the 90% C.L.U.L (statistical errors only) 
and an arrow pointing down. See Section [ill El for details. Dashed lines are the predictions for each sign from a data-driven 
hadron generator normalized to the measured points. The /dof for these fits are quoted. See Section Till PI for details. The 
width of the lines corresponds to an^ (see Table IVljl . the error on the ratio of free-decay muons to inclusive muon candidates. 
Black bands at the left edge of each panel show the pr-independent systematic error on the inclusive muon candidate yield. 
Shaded bands on each point show the systematic errors that affect the pr shape of the inclusive muon candidate spectrum. 
These last two systematic errors (Table [3 need to be included in the total error budget for the yield of free-decay muons, 
ono/Nb fTable lVIt . but are displayed separately since they are common to all components of the inclusive muon candidate 
yield, see Equations ^1 and ^1 Systematic errors are discussed in Section UlI II 
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FIG. 8: Points (measured data) show the longitudinal mo- 
mentum, measured at the vertex (p"*^), of particles that stop 
in MuID gap 3. The sharp peak is due to muons which range 
out in the absorber plate between gaps 3 and 4. The his- 
togram (Monte Carlo) shows the longitudinal momentum of 
all particles that stop in MuID gap 3 and do not decay before 
the pre-MuTR absorber. The Monte Carlo is normalized to 
the data for p"*'^ > 3GeV/c. Particles beyond the peak form 
a relatively pure sample of hadrons. 



Results of these fits are shown in Figure 

With these values for Xi{pT) and the hadron gener- 
ator input spectra, we could directly apply Equation 
to obtain the final punchtrough spectra. However, we 



made one further correction, described below, after find- 
ing that our assumption of exponential absorption does 
not hold when applied to GEANT simulations of the 
punchthrough process. 

Using our GEANT-based PHENIX simulation pro- 
gram, we generated data sets with both the FLUKA |^ 
and GHEISHA 69] hadronic interaction packages. In- 
put spectra for these data sets were given by our decay 
hadron generator. We selected all particles which did 
not decay before the pre-MuTR absorber. "Truth" val- 
ues for the punchthrough and gap 3 exclusive yields were 
obtained by splitting those particles based on the absence 
(gap 3 exclusive) or presence (punchthrough) of associ- 
ated charged particles with E4 > 100 MeV in MuID gap 
4. We varied E4 from 50—300 MeV and saw no significant 
change in the results. 

Using the known input spectra, known values for 
L3^4{0), and truth values for gap 3 exclusive yield, we ap- 
plied Equation [7| to the Monte Carlo data sets to obtain 
\iPT)- Due to statistical limitations we integrated our 
results over 9 and into two pT bins: 1 < pt < 2GeV/c 
and Pt > 2GeV/c. Values extracted for X^{pt) for the 
different hadronic interaction packages are listed in Ta- 
ble These values are consistent with those found for 
our measured data, listed above. 

Inserting these values for X^{pt) into Equation we 
obtained a prediction for Np{pT, 0)- Ratios of the truth 
values and predicted values for the punchthrough yield 
(i?7V^(py-)) for the different hadronic interaction packages 
are listed in Table Hvl One can see that these ratios devi- 
ate significantly from unity and that the two hadronic in- 
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FIG. 9: Points show the pr spectrum of (left) positively and (right) negatively charged hadrons stopping in MuID gap 3 ("gap 
3 exclusive yield") with statistical errors. Open diamonds show a power-law fit to the data, effectively a smoothed version 
of the data. Dashed lines are fits to the smoothed data using the hadron generator (normalized by the measured free-decay 
spectrum, as shown in FigureEJ and Equations |H1 and El to obtain values for the species-dependent nuclear interaction lengths, 

y{PT). 



TABLE III: Nuclear interaction lengths, A'(pt), for differ- 
ent particle species and pr bins (in GeV/c) for FLUKA and 
GHEISHA. Statistical errors on these values are « 2 mm. 



FLUKA GHEISHA 
A'(pt) [cm] A'(pt) [cm] 



Species 


1 < PT < 2 


PT > 2 


1 < PT < 2 


PT > 2 


n+ 


19.6 


24.5 


16.0 


21.1 


7r~ 


19.4 


24.8 


15.0 


19.3 


K+ 


24.4 


29.6 


24.9 


30.8 


R- 


20.5 


24.2 


17.2 


21.2 



teraction packages disagree on the direction of the devia- 
tion: the exponential absorption model tends to overpre- 
dict the punchthrough yield for FLUKA (i?^f ^^^^ < 1) 

and underpredict it for GHEISHA [R^V^^pJ]"'^ > 1). 
Relatively little data exists in the relevant momentum 



range that would allow us to conclude which, if either, 
of the hadronic interaction pac kage is correct. Measure- 
ments by RDIO and RD45 [T^ of the penetration depth 
of identified hadrons found that GHEISHA did well for 
protons and FLUKA did not. But neither did well for 
pions and no data exists for kaons. Furthermore, the re- 
sults are sensitive to the definition of a penetrating par- 
ticle: For RD 10/45 an incoming particle with any associ- 
ated charged particles in the 120 x 120 cm^ detector area 
at a particular depth was defined to have penetrated to 
that depth. In our measurement we reconstruct parti- 
cle trajectories and MuID hits are not associated with a 
road unless they are within a narrow search window sur- 
rounding the projected trajectory. Thus we are relatively 
insensitive to the leakage of a showering hadron. 

As a result of these uncertainties on the applicabil- 
ity of our exponential absorption model we incorporate 
a species and pT-dependent correction factor to Equa- 
tion 



N],{pT,e) = c\pT)K,,ipT, e) eM-Lmi\\pT)). (lo) 

The correction factors for pions and kaons are obtained from the average of the values of Rj^ji i^p^^ for the two 
packages, {RY{pt) = (^^^(^rf + ^n^M^^^^^ ^^^^^ ^'^^ ^^^^^'^ ™ Table Hvl The values of (RHpt) for a given 
species are not the same for the different pr bins, so we assume the values are valid at the average px of each bin 
{pt — 1.25GeV/c and 2.31GeV/c respectively) and use a linear extrapolation in px to obtain the final correction 
factors: 

C\pt) = {RTil <PT< 2GeV/c) + ((i?)^(pT > 2GeV/c) - (i?)*(l < pr < 2GeV/c))^^^[^|^^/^^^^ (11) 

I 

We assume that p's and p's have the same correction factors as the corresponding sign pions. Since p's and 
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TABLE IV: Ratios, Rn' (pt)^ truth values for the 
punchthrough hadron yield to those predicted assuming expo- 
nential hadron absorption for different particle species and pr 
bins (in GeV/c), for FLUKA and GHEISHA. Average values 
of the ratios for the two different hadronic interaction pack- 
ages, {R)^{pt), are smoothed across the pr bin at 2GeV/c 
to obtain correction factors for the exponential absorption 
model. Statistical errors on these quantities are ~ 10%. The 
maximum fractional difference in the ratios for the two differ- 
ent packages (32%) is incorporated into the systematic error 
estimate, as shown in Table IVni 



Species 


i < pt < ^ 


Pt > ^ 


Description 


TT 


U.lb 


U.oD 




TT 


0.91 


0.75 




K+ 


0.91 


1.00 


TjFLUKA 


R- 


1.17 


1.06 




7r+ 


1.48 


1.04 




7r~ 


1.47 


1.09 




K+ 


1.31 


1.07 


ryGHEISHA 


K' 


2.21 


1.69 




7r+ 


1.12 


0.95 




7r~ 


1.19 


0.92 




K+ 


1.11 


1.04 


{RYiPT) 


R- 


1.67 


1.38 




7r+ 


32% 


10% 




7r~ 


24% 


18% 




K+ 


18% 


3% 


SRm;m/C'{pt) 


K- 


32% 


22% 





p's make only small contributions to the punchthrough 
hadrons this simplifying assumption has little conse- 
quence. We incorporate the maximum fractional differ- 
ence in the ratios for the two packages (32%) into our 
systematic error estimate, as listed in Table IVTll 

We use Eauation llOl with the tabulated correction fac- 
tors, particle yields at the vertex given by our normalized 
hadron generator, the known value of Li{9), and the val- 
ues for A*(pt) determined from the measured gap 3 ex- 
clusive yield, to obtain the pr spectrum of punchthrough 
hadrons, Np{pt), as shown in Figure^] Systematic er- 
rors shown in this figure are discussed in Section llll H and 
listed in Table lyTTl We multiply Np{pt) by the fraction 
of the accepted z^tx range represented by each z^tx bin 
to finally obtain Np{pT, Zytx)- 

G. Background Tracks 

The main source of tracks which are not accounted 
for in the yield of punchthrough hadrons and free-decay 
muons, and which are not due to vertex- independent 
muons, are light hadrons which penetrate through the 
pre-MuTR absorber, decay into a muon, and are still re- 
constructed as a valid track. 




1 1.5 2 2.5 3 
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FIG. 10: Yield of positively (dotted line) and negatively (solid 
line) charged hadrons which penetrate to MuID gap 4. The 
curves are obtained from Equation 1101 as described above. 
Solid bands on the y-axis show the relative normalization er- 
ror on the inclusive muon candidate yield, a"^™ /Nj, see Ta- 
ble [V| This needs to be included in the total error budget 
for the yield of punchthrough hadrons, aNp/Np f Table IVTTl . 
but is displayed separately since it is common to all compo- 
nents of the inclusive muon candidate yield, see Eauations ll2l 
and 1131 Hatched bands on the j/-axis show cf^r™ (see Ta- 
ble lVIH , the normalization error on the ratio of punchthrough 
hadrons to inclusive muon candidates. The relative fraction 
of positive and negative punchthrough hadrons can move up 
and down together by this amount. Shaded bands around the 
extracted punchthrough hadron yield show the systematic er- 
rors on cr^^ which can affect the pr shape of the relative frac- 
tion of positive and negative punchthrough hadrons (positives 
and negatives can move independently) . These are dominated 
by differences in the applicability of the simple exponential 
absorption model observed for FLUKA and GHEISHA. Sys- 
tematic errors are listed in Tables [V| and IVIII and discussed 
in Section ITTTTl 



A simulation of single pions thrown into the muon 
arm acceptance shows that the number of hadrons which 
decay after the pre-MuTR absorber and penetrate to 
MuID gap 4 is only 5-10% (increasing with increasing 
Pt) of the z„fj;-averaged number of free-decay muons, 
NdIpti Zytx = 0). This ratio will be suppressed by the 
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fact that tracks which decay are less likely to be recon- 
structed successfully. It is further suppressed by our 
punchthrough calculation procedure: the number of such 
tracks which stop in MuID gap 3 is roughly half the num- 
ber that penetrate to gap 4; these will be counted in our 
calculation of the punchthrough hadron yield. 

We express our estimate for the yield of background 
tracks not otherwise accounted for as Nb{pt) = 5% x 
Nd{pt, Zvtx = 0). The systematic uncertainty assigned 
to this quantity, ±5% x Nd{pt, Zvtx = 0), covers the 
extreme possibilities that the Nb/Njj is unsuppressed 
or fully suppressed by reconstruction and punchthrough 
procedures, as described above. 

This estimate was verified in a simulation of tt^'s and 
K~^s which were thrown into the muon arm acceptance 
and fully reconstructed. The reconstructed track infor- 
mation, together with the Monte Carlo truth informa- 
tion, allows us to eliminate uncertainties due to mis- 
reconstruction of the track pT and due to determination 
of whether a track which penetrated to the last gap did 
so in a reconstructible fashion. Systematic errors on this 
estimate are discussed in Section UlI II 



H. Vertex-Independent Muons 

Figure ^2 shows the yield of inclusive muon candi- 
dates, Ni(j>T, Zytx), with contributions from individual 
components (free-decay muons, punchthrough hadrons, 
and background tracks) shown as well as their sum. The 
vertex-independent muons can be seen as the clear excess 
above the calculated background sources. The systematic 
error bands shown on the component sums are discussed 
in Section inni and hsted in Tables IVll and Ivlll 

We obtain the yield of vertex-independent muons by 
applying Equation |21 in each px bin, subtracting the 
hadronic contributions from the inclusive muon candi- 
date yield, and averaging over z^tx bins. This yield is 
shown, before averaging over Zytx to demonstrate the ex- 
pected vertex independence, in Figure [T^ 

We make one final correction for momentum scale. The 
observed mass of the J ftp, reconstructed with the same 
code and in the same data set, is higher than the nominal 
value by « 100 MeV(3%')[5ll|. However, in a higher statis- 
tics data set the momentum scale accuracy is verified to 
within 1% by our observation of the accepted value for 
the mass of the J/i/' Also, the peak observed in the 
longitudinal momentum distribution of particles stopping 
in MuID gap 3 (see Fig.|SJ) is within 0.5% of the predicted 
value. We therefore assume that the momentum scale is 
high by 1.5% (splitting the difference between and 3%). 
This results in a momentum scale correction factor to the 
prompt muon yield of 0.94 + 0.987 x (pT[GeV/c] - 1). We 
assume a 100% systematic error on this correction factor, 
as shown in Tab. 

Finite momentum resolution can cause a similar ef- 
fect. Contributions from energy loss fluctuations, mul- 
tiple scattering and chamber resolution combine to give 



5p/p « 5% for the momentum range used in this anal- 
ysis. Finite resolution, combined with an exponentially 
falling spectrum, artificially hardens the measured spec- 
trum. For 1 < Pt < 3GeV/c this hardening increases 
the normalization of the yield by « 3.7%. However, this 
is accounted for in our efficiency determination since we 
use a realistic pt spectrum as input. Therefore we apply 
no explicit correction for this effect. 

The final values for the vertex-independent muon cross 
section, obtained from Equation |31 are shown in Fig- 
ure 1131 Points in this figure have been placed at the 
average px value of the bin contents to account for bin 
shifting in the steeply falling distributions. Systematic 
errors shown in this figure are discussed in Section IIII II 
and hsted in Tables Ivl - IVlTl 

I. Systematic Errors 

Many sources of systematic error on the yield of vertex- 
independent muons, iV^, are common to the different 
components of the inclusive muon candidate yield. In 
order to account for this we rewrite Equation |21 (making 
the Pt and Zytx dependencies implicit) as: 

^ NiX {N^/Ni) (12) 
= Nix {{Ni -Nd-Np- Nb)/Ni) 
= Ni X (1 - Rd - Rp - i?s), 

where Rj — Nj/Nj is the fraction of the inclusive muon 
candidate yield attributed to the j*'* component. We can 
now write the systematic error on N^^ as: 

^iv. - ^{<yNr/NjrN^, + (4, + 4p + (13) 

(TAT^, as determined below, is displayed in Figures 1111 
andini Note that the errors for positives are significantly 
larger than for negatives. This is due to the much larger 
relative contribution to positive inclusive muon candi- 
dates from punchthrough hadrons, which is due to the 
relatively small size of the K'^ nuclear interaction cross 
section. 

Error sources contributing to crjVj are quantified in Ta- 
ble Error sources contributing to aNo and ap^ are 
quantified in Table lVll Error sources contributing to ajsip 
and (Tpp are quantified in Table IVlTl Note that in these 
tables we separately list errors that affect the overall nor- 
malization (cr/iV"""") and the shape of the pr spectrum 
(a/NP-^ ). The error on Rb is taken to be 100% of its esti- 
mated value: ap^ = Nb/Ni = 0.05 x Noizytx = 0)/iV/. 

Values for aNj/Nj are displayed in Figure Values 
for ap^ and app are displayed in Figures ^1 and [7| re- 
spectively. We insert cfnj/Nj, ap^, ctrp, and apg into 
Equation 1131 as part of the final systematic error on N^. 

To get the vertex-independent muon cross section, as 
defined in Equation|3|and displayed in Figures lT^ and lTBl 
we need to add in quadrature the errors on A^^, o^bbc^ 
s'bbc^ ■ The error on iV^ is obtained from the components 
above according to Equation 1131 As mentioned 
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FIG. 11: Points show the yield of (left) positively and (right) negatively charged inclusive muon candidates vs. z^tx for different 
Pt bins with statistical errors. Dotted, solid and dashed lines show contributions from decay muons, punchthrough hadrons 
and background tracks, respectively. Shaded bands show the systematic error around the sum of these components, as listed 
in Tables rVlVni and discussed in Section [111 II 
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FIG. 12: Points show the yield of (left) positively and (right) negatively charged vertex-independent muons vs. Zvtx for different 
Pt bins with statistical errors. The dashed lines show the yield for each pt bin, averaged over z^tx- The shaded bands around 
show the systematic error on the sum of the contributions to the inclusive muon candidate yield from light-hadronic sources, 
as listed in Tables IVl Vlll and discussed in Section fill II 
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FIG. 13: pt spectrum of vertex-independent muons. Error 
bars indicate statistical errors. One point with unphysical 
(less than zero) extracted yield is shown as an arrow pointing 
down from the 90% C.L.U.L. Shaded bands indicate system- 
atic errors, as listed in Tables IVIVllI and discussed in Sec- 
tion Eni 



above, c^'^p was determined by a van der Meer scan to 
be 21.8 mb, with an error of 9.6% . We assign an error 
of 5% to e^BcT' which was found to be 0.75 through two 
different methods IsHIBOI. 



IV. CHARM CROSS SECTION 

The charm production cross section obtained from the 
yield of vertex-independent muons (or from the yield 
of non-photonic electrons, or D mesons) is necessarily 
model dependent since we do not measure the charm 
quarks directly. We use PYTHIA to convert our mea- 
surement of the vertex-independent muon yield into an 
estimate of the differential charm production cross sec- 
tion at forward rapidity, dacc/dy\y=i.G-, in a procedure 
very similar to that used in PHENIX measurements of 
charm production at?; = o[l|i|i,[ll[ll[ll[l3. We 
use PYTHIA version 6.205 with parameters tuned to re- 
produce charm production data at SPS and FNAL |23| 
and single electron data at t he ISR [llllllzi- Tuned 
parameters are listed in Table IVTllI The meaning of each 
parameter is more thoroughly defined in the PYTHIA 
manual [TSf. 

Vertex-independent muon sources, predicted by a 
PYTHIA simulation using the same parameters (except 
that MSEL is set to 2 to generate unbiased collisions), 
are listed in Table IIXI These sources include decays 
of hadrons containing a heavy quark, and light vec- 
tor mesons with a decay length too short to be mea- 
sured with the existing experimental apparatus (p, lu, (f). 
Their pT spectra are shown in Figure 1141 Contributions 
from quarkonium decays, Drell-Yan and r lepton de- 
cays are negligible. This shows that vertex-independent 
muon production in our acceptance is dominated by 



muons from decay of charm hadrons, although for pT > 
2.5GeV/c the contribution from decays of hadrons con- 
taining a bottom quark is starting to become important. 

charm. Bottom, and Light Vector Meson Contributions in PYTHIA 
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FIG. 14: PYTHIA calculation showing the major contribu- 
tions to the vertex-independent muon pr spectrum. Solid, 
dashed and dotted lines show the yield from charm, bottom 
and short-decay length light vector mesons {p,u},(j>) respec- 
tively. 

This simulation also gives the distribution of charm 
quarks (pt vs. y) that produce a muon in our accep- 
tance, as shown in Figure 1151 This demonstrates that 
the vertex-independent muons we measure sample charm 
quarks down to px ~ IGeV/c, over a narrow rapidity 
slice centered at y — 1.6. 
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FIG. 15: PYTHIA results for the pr vs. y distribution (linear 
z-scale) of charm quarks that produce a muon in the PHENIX 
acceptance. 

Figure [TCI shows a comparison of the measured vertex- 
independent negative muon spectrum (from Figure I13() 
to the prediction of this default PYTHIA simulation and 
to a FONLL calculation [23, ll^l • One can see that the 
measured values significantly exceed both predictions. 
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TABLE V: Sources of systematic error on the calculation of Ni, the yield of inclusive muon candidates. 



Error Source cr/A^" 



Comment 



Momentum Scale 6.0% 

£acc 10% 

Erec 9.0% 

£u.ar 5.0% 

etrig 4.7% 

aNj/Ni 16.3% 



(pT[GeV/c] - 1) X 1.3% Taken to be 100% of the correction. 

(pT[GeV/c] - 1) X 1.5% Taken to be 20% of the correction. 

Taken to be 25% of the correction. 

(pT[GeV/c] - 1) X 5.0% Error on ct"°"" taken to be 20% of the correction. 

Error on a'"^ taken to be maximum of observed pr variation. 

Taken to be the difference in etrig obtained with different procedures. 

(pT[GeV/c] — 1) X 5.4% Add all rows in quadrature. 



TABLE VL Sources of systematic error on Rd, the ratio of free-decay muons to inclusive muon candidates, and Nd, the 
absolute yield of free-decay muons. 



Error Source 



a/N" 



Comment 



Decay flight path 5% 

Zvtx fit range 3.3% 
Input hadron spectrum 

Decay normalization 7% 
9.2% 

aNjNi 16.3% 

aNo/No 18.7% 



Variation due to shift of \d by ±3 cm. 

Variation seen using \zvtx\ < 20cm, \zvtx\ < 40cm. 

(Pt [GeV/c] — 1) X 5.0% Spectrum normalized to data, so the only 
uncertainty here is in the pr dependence. 

Statistical uncertainty in fit to observed free-decay muon yield. 

(pT[GeV/c] — 1) X 5.0% Add all previous rows in quadrature. 

(pT [GeV/c] - 1) X 5.4% From Table |V| 

(pt [GeV/c] — 1) X 7.4% Add previous two rows in quadrature. 



TABLE VII; Sources of systematic error on Rp, the ratio of punchthrough hadrons to inclusive muon candidates, and Np, the 
absolute yield of punchthrough hadrons. 



Error Source 



Comment 



Exponential absorption model 



p3 

'-scale 



p and p contributions 



N^^°P normalization 



(TNr/Ni 
(^Np/Np 



23% 

10% 

10% 
27% 
16.3% 
31.5% 



32% 











Maximum fractional difference between C^{pt) 
for FLUKA and GHEISHA, from Table HVl 
Dominated by ambiguity in the definition of which 
particles should be reconstructed. 
Variation with extreme assumptions on the p and 
p nuclear interaction length. 

Statistical uncertainty in fit to observed gap 3 exclusive yield. 
Add all previous rows in quadrature. 
(pT[GeV/c] - 1) X 5.4% From Table|V| 

~ 32% Add previous two rows in quadrature. 




32% 



The spectrum also appears to be somewhat harder than 
the PYTHIA spectrum vi/ith the parameters hsted in Ta- 
ble EIIIl 

We scale the charm (only) contribution to the 
PYTHIA vertex-independent muon px spectrum such 
that the total spectrum (including the small contribu- 
tions from open bottom and vector mesons) matches the 
central values of the measured vertex-independent nega- 
tive muon spectrum. Only statistical errors are used in 
the fit. Note, larger systematic errors for the positive 
muon spectrum preclude a significant measurement for 
that charge sign. We multiply the scale factor from the 
fit (2.27) by the PYTHIA value for the charm produc- 
tion cross section, dace/ dy\yXi Q^^ (O.lOTmb), to obtain 
dacc/dy\'y^^^i^'^ = 0.243 ± 0.013(stat.) mb. 



We distinguish between two different sources of sys- 
tematic uncertainty on the extraction of the charm cross 
section: 1) uncertainty in the PYTHIA calculation, and 
2) uncertainty in the data, which is largely independent 
of PYTHIA. 

We determined the uncertainty in the data (±43%) 
by refitting PYTHIA to the data at the minimum and 
maximum of the la systematic error band. 

We determined the uncertainty in the PYTHIA calcu- 
lation with a systematic study in which we varied sim- 
ulation parameters, extracted the new simulated vertex- 
independent negative muon spectrum, normalized to the 
measured spectrum, and extracted dacc/dy\y=i^i^^'^ for 
the modified parameter sets. We varied PDF libraries, 
the hard scattering scale, the charm quark mass, the 
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TABLE VIII: Tuned PYTHIA parameters (default settings 
for this analysis) for determination of charm production cross 
section central value. 

Parameter Value Meaning 



p+p ^ii + X(y = -1.6.-^ = 200 GeV ) 



MSEL 

MSTP(32) 
MSTP(33) 
MSTP(52) 
MSTP(51) 
MSTP(91) 
PARP(31) 
PARP(91) 
PARP(93) 
PMAS(4,1) 
D+/D° 



4 Heavy quark production every event 

(gluon fusion + q/q annihilation). 

4 Hard scattering scale, = s. 

1 Use ii'-factor. 

2 Use PDF libraries. 

4046 Select CTEQ5L PDF hbraries [76]. 

1 Use Gaussian distribution for intrinsic kr- 

3.5 if-factor. 

1.5 <kT> (GeV/c). 

5.0 Maximum kr (GeV/c). 

1.25 nic (GeV/c). 

0.32 Default charm chemistry ratio. 



TABLE IX: Percentage contribution of different sources of 
vertex-independent muons within our acceptance (1 < pr < 
3GeV/c and 1.5 < \ti\ < 1.8), from PYTHIA, with parame- 
ters listed in Table IVIIII (except that MSEL = 2 to generate 
minimum bias collisions). 



Source 



Contribution 



Open charm 
Open bottom 

p,UJ,<j) 

Quarkonia 
Drell-Yan 
r leptons 



84.6% 

6.9% 

8.1% 

< 0.1% 

< 0.1% 
0.4% 



intrinsic /ct value, the D+/D° ratio, charm production 
mechanism selections, and open bottom and vector me- 
son scaling assumptions. The parameter sets used and 
the results of this study are summarized in Table Ixl 

The PYTHIA charm cross section varies substantially 
{A{dacc/dy)\yXj^i^^'^) « 4) for the chosen parameter 
sets. However, the extracted experimental charm cross 
section is relatively stable {A{daac/dy)\i^J^i%^^^ < 0.36). 
This is due to the fact that the parameter set changes 
have relatively minor effects on the shape of the pre- 
dicted vertex-independent muon pT spectrum, and we 
obtain the experimental charm cross section by normal- 
izing the PYTHIA charm cross section by the ratio of the 
measured and predicted muon pT spectra. 

One way to visualize this is to plot (see Figure I17() 
the vertex-independent muon yield in our acceptance 
per event in which a cc pair is created for the differ- 
ent PYTHIA parameter sets. Due to our procedure, pa- 
rameter sets which give similar vertex-independent muon 
yields per cc event in the low pT region (which dom- 
inates the fit) will necessarily give similar values for 
dacc/dy\y^i%^'^ , whatever the PYTHIA charm cross 



K lo- 



ll 

O 10 

I 





2.5 3 
p (GeV/c) 



FIG. 16: The top panel shows the measured pr spec- 
trum of vertex-independent negative muons from Figure [T^ 
the PYTHIA prediction using settings listed in Table IVIIII 
without scaling the charm contribution (dotted line), and 
a FONLL calculation (solid line with systematic error 
band) ^ol . Fz^ . The bottom panel shows the ratio of the 
measured spectrum to the FONLL calculation with statistical 
(error bars) and systematic (bands) uncertainties on the data, 
as well as the theoretical uncertainty (shaded band around 1). 
The dashed hue shows the PYTHIA/FONLL ratio. 



section is. 

The largest variation in the predicted muon yield 
at pt = IGeV/c per cc event is seen for simula- 
tions in which the intrinsic kx is varied from its de- 
fault value (< kr >= 1.5GeV/c) to the value ex- 
pected from arguments based on Fermi momentum (case 
4c, < kT >= 0.3GeV/c), or to a value which best 
reproduces the measured spectrum at higher px (case 
4d, < kx >= 3.0GeV/c). These parameter sets also 
result in the largest variation in dace/ dy\y^^^^-^ , as 
shown in Table |X] We use the cross section values ob- 
tained in this pair of simulations to define the system- 
atic uncertainty in our measurement due to the uncer- 
tainty in our PYTHIA calculation. This gives us our 
final answer: dace/ dy\y^i,Q = 0.243 ± 0.013(stat.) ± 



0.105(data syst.) tj^;^^? (PYTHIA syst.)mb. 

Figure ^| shows the PHENIX charm rapidity spec- 
trum. The result of this analysis (mirrored about y = 
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FIG. 17: (Color online) PYTHIA results with different parameter sets for the negative vertex-independent muon pr spectrum 
per event in which a cc pair is created. The solid line in each panel shows the result when using default settings listed in 
Table IVrm Legends indicate the correspondance between line style and the simulation case label. Parameter sets for each case 
label are given in Table DO 



since this is a symmetric collision system) is plotted along 
with the result for dace/ dy\y^Q uM- order to compare 
with the data at y = the systematic uncertainty on the 
data from this analysis is shown as the quadrature sum of 
the two sources of systematic uncertainty described above 
(data and PYTHIA). Theoretical curves from PYTHIA 
(case 1 and case 5a)^0NLL jMEl, and an NLO cal- 
culation from Vogt [831 are also displayed. 

In the top panel of the figure PYTHIA with the de- 
fault parameter set (Case 1) is fit to the two PHENIX 
points with statistical and systematic errors added in 
quadrature. Other theory curves are normalized so 
that they are equal at y — in order to allow 
shape comparisons. As shown in Table H different 
PYTHIA parameter sets differ in the predicted ratio 
da,,/dy\PIli''^/dcj,-Jdy\P^™^ by > 30%. Unfortu- 
nately, current systematic error bars preclude any con- 
clusions about the charm production rapidity shape. 

In the bottom panel of the figure the theory curves 
are unnormalized to allow an absolute comparison. The 
quoted theoretical uncertainty bands for the FONLL 
and NLO calculations are also shown. We note that, 
although our data are above the FONLL prediction, 
the error bars touch. This is in contrast to the sit- 
uation for the vertex-independent muon cross section, 
shown in Figure 1161 where the data are significantly 



above the prediction. The larger disagreement in the 
vertex-independent muon cross section is presumably due 
to different treatment of the fragmentation process in 
PYTHIA and FONLL HSIlllll. 

V. CONCLUSION AND OUTLOOK 

We have measured muon production at forward ra- 
pidity (1.5 < I77I < 1.8), in the range 1 < pr < 
3GeV/c, in = 200 GeV p + p collisions at RHIC. 
We determined and subtracted the contribution from 
light hadron sources {tt,K,p) to obtain the vertex- 
independent muon yield which, for the pt range mea- 
sured in this analysis, and in the absence of new 
physics, arises dominantly from the decay of D mesons. 
We normalized the PYTHIA muon spectrum resulting 
from the production of charm quarks to obtain the 
differential cross section for charm production at for- 
ward rapidity: dacc/dy\y=i.6 — 0.243 ± 0.013(stat.) ± 
0.105(data syst.) tf5:^87 (PYTHIA syst.) mb. This is com- 
patible with PHENIX charm measurement at y = 0, al- 
though even further above predictions from PYTHIA and 
FONLL. Large systematic uncertainties in the current 
measurement preclude statements about the rapidity de- 
pendence of the charm cross section. 
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TABLE X: Results for PYTHIA simulations with different parameter sets used to explore the systematic error on the charm 
cross section due to model uncertainties. The top of the table details the different parameter sets tested. Unless otherwise 
noted, parameters are the same as those listed in Table fvTTTl The bottom of the table gives the results for different simulations: 
The 1"* column identifies the simulation; the 2"'' column gives the total charm production cross section given the chosen 
PYTHIA parameter set; the 3'^'' column gives the differential charm production cross section at j/ = 1.6; the 4*'' column gives 
the normalization factor needed to fit the PHENIX data; the S**^ column gives the differential charm production cross section 
at y = 1.6 for PHENIX data (the product of the 3'^'' and 4*'' columns); the 6*'' column gives the fractional difference between 
the results for each simulation compared to the simulation with the default PYTHIA parameter set; the last column gives the 

4-- J /J iPYTHIA II I, iPYTHIA 

ratio dacc/dy\y^i,a /dacc/dy\y^a 



Case PYTHIA Settings 

1 Default settings, see Table IVTm 



2a 


MSTP(51) = 


: 4032, CTEQ4L PDF hbraries [78]. 








2b 


MSTP(51) = 


: 5005, GRV94LO PDF libraries [79]. 








2c 


MSTP(51) = 


: 5012, GRV98LO PDF libraries [80]. 








2d 


MSTP(51) = 


: 3072, MRST (c-g) PDF libraries [81j. 








3a 


MSTP(32) = 


: 1, Q2 = 2siu/{f + P + u^). 








3b 


MSTP(32) = 


--2,Q^ =p?r + (mi + ml)/2. 








3c 


MSTP(32) = 


: 3, = min(— f, —u). 








3d 


MSTP(32) = 


= 5, q2 = -i. 








4a 


PMAS(4, 1) : 


= mc = 1.15GeV/c. 








4b 


PMAS(4, 1) : 


= mc = 1.35 GeV/c. 








4c 


PARP(91) = 


< kr >= 0.3GeV/c. 








4d 


PARP(91) = 


< kT 3.0GeV/c. 








4e 


MSTP(68) = 
PARP(67) = 


: 2, Maximum virtuality scale and matrix element matching 
4, Multiplicative factor applied to hard scattering scale. 


; scheme. 




5a 


PARP(31) = 


if- fact or = 1, 










MSEL = 1, Hard scattering enabled. 








5b 


PARP(31) = 


if- fact or = 1, 










MSEL = 1, Hard scattering enabled. 










All other parameters untuned. 








6 


D+/D° = 0.45 [82j. 








7 


Open bottom and vector mesons scale with charm. 








Case 


PYTHIA 


dac.c./dy\yZ-i.6'^ Normalization 


da.,ldy\'^S,T''' 


^da^./dyl^yS.T"' 


da.,/dyfyZ™^/ 




(mb) 


(mb) to Data 


(mb) 


(%) 




1 


0.658 


0.107 2.27 


0.243 


- 


0.67 


2a 


0.691 


0.111 2.10 


0.232 


-4.5 


0.69 


2b 


0.698 


0.112 2.09 


0.233 


-3.9 


0.71 


2c 


0.669 


0.109 2.18 


0.238 


-1.7 


0.73 


2d 


0.551 


0.088 2.67 


0.236 


-2.9 


0.71 


3a 


1.520 


0.243 1.12 


0.271 


11.8 


0.84 


3b 


0.863 


0.139 1.63 


0.226 


-6.7 


0.71 


3c 


1.501 


0.242 1.11 


0.267 


10.2 


0.84 


3d 


1.104 


0.178 1.45 


0.258 


6.4 


0.78 


4a 


0.905 


0.145 1.73 


0.252 


3.7 


0.67 


4b 


0.487 


0.078 2.91 


0.226 


-6.7 


0.64 


4c 


0.658 


0.104 2.81 


0.292 


20.4 


0.66 


4d 


0.658 


0.104 1.50 


0.156 


-35.8 


0.63 


4e 


0.658 


0.106 2.09 


0.220 


-9.2 


0.63 


5a 


0.435 


0.068 3.91 


0.266 


9.4 


0.80 


5b 


0.385 


0.058 4.67 


0.271 


11.7 


0.79 


6 


0.658 


0.107 2.38 


0.255 


5.0 


0.67 


7 


0.658 


0.107 2.20 


0.236 


-2.9 


0.67 
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PYTHIA ■ 
(Case 5a) 
FONLL ■ 
NLO 



tics, and with hadrons stopping in MuID gap 2, will allow 
a completely data-driven approach to the calculation of 
the punchthrough yield. This will eliminate the reliance 
on hadronic interaction simulation packages, differences 
in which are the largest source of systematic error at 
high pt- Analogous measurements are also being carried 
out for d + Au, Cu + Cu, and Au + Au 84] collisions 
at y^SNN = 200 GeV. These will allow determination of 
the magnitude of nuclear modification effects on charm 
production at forward rapidity. 
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PYTHIA - 
(Case 1} 

PYTHIA ■ 
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NLO E 




FIG. 18: (Color online) Comparisons of measured charm ra- 
pidity distributions, docc/dy vs. y, to theoretical predictions. 
Data points at y = ±1.6 are from this analysis (the point at 
1/ = 1.6 is reflected through y — 0). The point at y = is the 
PHENIX measurement of charm through semileptonic decay 
to electrons Error bars on the data points indicate sta- 
tistical uncertainties and boxes indicate systematic uncertain- 
ties. The top panel shows rapidity spectra from two PYTHIA 
parameter sets (see Table^for details), FONLL 113,173, and 
an NLO calculation [gl. The PYTHIA curve with the de- 
fault parameter set (Case 1) was fit to the two PHENIX data 
points with statistical and systematic errors added in quadra- 
ture. All other theory curves were normalized so that they 
are equal at j/ = to allow shape comparisons. The bottom 
panel shows the theory curves unnormalized. Theoretical un- 
certainties associated with the FONLL and NLO calculations 
are indicated with shaded bands. 



The systematic uncertainty in the data is dominated by 
uncertainty on the determination of the fractional con- 
tribution of decay muons. This will be improved with 
higher statistics data sets (already collected) which will 
allow better measurements of the Zytx dependence of par- 
ticle production. Final results for identified particle pt 
distributions in p + p collisions by BRAHMS will also be 
invaluable for improving the input to our hadron gener- 
ator. The systematic uncertainty in PYTHIA is domi- 
nated by differences observed when the intrinsic < fcy > 
is varied. In order to reduce this uncertainty we need 
to reduce the allowed parameter space by improving 
the measurement of the high pT portion of the vertex- 
independent muon spectrum, where the error is domi- 
nated by the uncertainty in the yield of punchthrough 
hadrons. Data sets (already collected) with higher statis- 
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